The  Effects  of  Near- Wall  Control  on  the  Structure  and  Dynamics 
of  Wall-Bounded  Turbulent  Flows 


Kenneth  Breuer 
Division  of  Engineering 
Brown  University 
Providence  RI  02912 


20060302  164 


This  report  is  the  final  technical  report  for  the  ONR  award  N00014-04- 1-0604  (Brown  Account  5-25030) 
Any  technical  questions  should  be  addressed  to  kbreuer@brown.edu 


Contents 


Table  of  Contents  ii 

List  of  Tables  iii 

List  of  Figures  iv 

1  Introduction  1 

1.1  Background .  1 

1.2  Review  of  scaling  arguments .  1 

1.3  Current  Experiments .  2 

2  Experimental  Setup  and  Procedure  3 

2.1  Actuator  Design  and  Fabrication .  3 

2.2  Flow  Facility .  6 

2.3  Particle  Image  Velocimetry  System .  7 

2.4  Experimental  Procedure .  7 

3  Results  and  Discussion  10 

3.1  Mean  Flow .  10 

3.2  First  Order  Statistics .  10 

3.2.1  Actuation  Data .  10 

3.2.2  Integrated  Flow  Rate  and  Production .  13 

3.2.3  Phase  Synchronization .  14 

3.3  Second  Order  Statistics  .  16 

3.3.1  Two-point  correlations .  16 

3.3.2  Conditional  Averaging .  17 

4  Conclusions  19 

Bibliography  21 


ii 


List  of  Tables 


2.1  Forcing  frequencies,  amplitudes,  dimensionless  numbers  and  flow  parameters  for  the  four  flow 

rates  used .  6 


iii 


List  of  Figures 


2.1  Schematic  of  the  Lorentz  Force  Actuator  showing  side  and  top  views  respectively.  The  Lorentz 
force  acts  in  the  spanwise  direction,  out  of  the  page  on  the  side  view  and  from  top  to  bottom 

on  the  top  view  [1] .  4 

2.2  Composite  photograph  of  Lorentz-force  actuator  showing  the  sub-structure  (flux  plate,  align¬ 

ment  pins  and  magnets)  and  the  electrode  layer  [1],  The  flow  in  this  picture  is  from  left  to 
right,  and  the  Lorentz  force  is  from  top  to  bottom.  The  streamwise  pitch  of  the  magnets  and 
electrodes  is  3.125mm .  4 

2.3  Photographs  showing  (a)  PCB  in  aluminum  mold  with  tape  border  and  1/8”  spacer  inserted 

into  corner  and  (b)  finished  plate  with  epoxy  coating  and  “levelling  feet”  embedded  on  the 
back  surface.  Holes  drilled  through  the  flux  plate  to  allow  seepage  are  also  visible .  4 

2.4  Schematic  of  flow  facility,  showing  the  back  pressure  and  water  supply  tanks,  development 

section,  test  area  with  actuator  and  return  flow  path .  5 

2.5  Schematic  showing  side  and  top  views  of  PIV  setup.  A  laser  sheet  in  the  x-y  plane  is  produced 

by  an  Nd:YAG  laser  mounted  above  the  channel.  The  image  size  is  approximately  9  x  11mm.  6 

2.6  Experimentally  measured  voltage  divided  by  current  plotted  with  actuation  frequency  ( T+ ). 

V/I  remains  approximately  constant  for  all  experiments .  7 

2.7  PIV  baseline  channel  velocity  profiles  showing  both  mean  flow  and  fluctuations.  Profiles 
demonstrate  fully  developed  turbulent  flow  at  Re*  =  167.  All  quantities  are  normalized  by 

the  friction  velocity,  uT .  8 

2.8  Velocity  profiles  at  a  higher  Reynolds  number  (Re*  =  330)  showing  both  the  mean  flow  as 

well  as  the  fluctuating  components .  8 

3.1  Baseline  unactuated  mean  profiles  plotted  with  numerical  turbulence  statistics  for  a  range  of 

Reynolds  numbers  from  Moser  et  al.  [2] .  11 

3.2  Mean  flow  profiles  for  a  fixed  forcing  amplitude  with  changing  frequency  at  Re*  =  167.  The 

fixed  amplitude  that  was  chosen  was  the  highest  possible  given  the  actuator  design.  Actuation 
both  increases  and  decreases  the  fluctuations,  with  a  local  minimum  around  T+ St/Re*  =  2.9.  11 

3.3  Changes  in  Reynolds  stress  and  turbulence  production  for  a  constant  amplitude  and  changing 

frequencies.  The  lowest  production  accompanies  the  highest  local  flow  rate  as  well  as  the 
greatest  reduction  in  shear  stress .  12 

3.4  Effects  of  forcing  on  the  mean  flow  are  shown  for  a  higher  flow  rate  (Re*  =  330).  One  forcing 

frequency  and  three  amplitudes  are  plotted  withthe  baseline .  12 

3.5  Changes  in  Reynolds  stress  and  turbulence  production  at  Re*  —  330.  The  turbulence  produc¬ 
tion  is  greatly  suppressed  around  the  ‘sweet  spot’  where  T+  St/ Re*  =  2.2 .  13 

3.6  Percent  change  in  flow  rate  and  production  terms  over  a  wide  range  of  forcing  frequencies  and 
amplitudes  for  a  flow  rate  of  Re*  =  167.  The  maximum  increase  in  flow  rate  is  accompanied 

by  the  greatest  decrease  in  turbulence  production,  near  the  ‘sweet  spot’ .  13 

3.7  Percent  change  in  flow  rate  and  production  terms  over  a  wide  range  of  forcing  frequencies  and 
amplitudes  for  a  flow  rate  of  Re*  =  334.  The  data  shown  represents  five  forcing  frequencies 

and  five  forcing  amplitudes .  14 


IV 


3.8  Changes  in  the  turbulence  production  and  flow  rate  as  functions  of  the  scaling  parameter, 

StT+/Re*.  Data  from  four  different  Reynolds  numbers  are  shown  and  a  wide  range  of  forcing 
amplitudes  and  frequencies.  Data  collapses  well,  showing  a  minimum  in  production  and 
maximum  in  flow  rate  at  T+St/Re*  between  2  and  3 .  15 

3.9  Mean  and  Reynolds  Shear  Stress  taken  at  four  different  phases  of  the  forcing  cycle  for  a 

Re*  =  200 .  15 

3.10  Turbulence  production  calculated  for  different  phases  of  the  forcing  cycle .  16 

3.11  RuU{yi,yi)  and  Rvv(yi,  2/2)  at  the  location  of  maximum  turbulence  production,  y+  =  12.  In  all 
actuated  cases,  the  correlation  decreases,  with  a  sweet  spot  at  approximately  StT+  /Re*  =  2.9  16 

3.12  Ruu(yuy2)  in  the  log  layer  for  (a)  Re*  =  167  at  y+  =  76  and  (b)  Re*  =  330  at  y+  =  122.  .  .  17 

3.13  Probabilities  plotted  as  a  function  of  the  threshold  parameter,  k,  for  a  frequency  sweep  at  a 

constant  forcing  amplitude  for  both  low  and  high  Reynolds  number  flow .  18 


v 


Chapter  1 


Introduction 


1.1  Background 

Improved  understanding  of  unsteady  flow  physics  has  led  to  a  wealth  of  flow  control  applications  in  which 
fluid  perturbations  are  deliberately  introduced  into  a  large-scale  flow  with  the  goal  of  affecting  its  global 
characteristics.  Research  into  these  methods  has  included  separation  control,  enhanced  mixing,  and  reduction 
of  turbulent  skin  friction.  Interest  in  turbulence  control  has  grown  due  to  numerical  computations  at  low 
Reynolds  number  flow  [3,  4,  5]  which  indicate  that  the  turbulence  and  drag  can  be  reduced  by  a  variety 
of  fluidic  actuations.  The  difficulty  with  these  methods  is  that  even  though  computationally  they  provide 
promising  results,  implementing  them  in  a  physical  experiment  proves  very  difficult.  Other  actuation  methods 
that  have  been  studied  experimentally  include  localized  heating  [6],  piezo-electric  flaps  [7],  oscillatory  blowing 
[8],  synthetic  jets  [9],  surface  motion  [10]  and  plasma  discharge  [11]. 

Given  the  previous  problems  with  designing  robust  actuators,  an  appealing  approach  for  turbulence 
suppression  in  an  electrically  conducting  fluid  is  Lorentz  force  actuation.  Lorentz  force  actuators  have  no 
moving  parts  and  thus  in  principle  are  easy  to  implement.  When  a  current,  I,  and  a  magnetic  field,  B, 
interact,  the  Lorentz  body  force,  Fl,  is  generated: 

Fl  =  IB  a  (1.1) 

where  a  is  the  characteristic  length  over  which  the  current  flows.  The  motivation  for  this  application  comes 
from  the  recognition  that  an  oscillatory  wall  motion  disrupts  the  spanwise  structure  of  the  equilibrium 
turbulent  flow,  disturbing  the  turbulence  production  cycle  and  consequently  the  turbulent  drag  [5,  12]. 

Recent  experimental  results  have  been  published  by  Di  Cicca  et  al.  [13]  and  Pang  and  Choi  [14]  document¬ 
ing  reductions  in  shear  at  the  wall  using  oscillatory  wall  motion  and  Lorentz  force  actuators,  respectively. 
Di  Cicca  et  al.  documented  particle  image  velocimetry  measurements  of  a  turbulent  boundary  layer  with 
an  oscillating  wall.  Experiments  were  conducted  in  a  closed  loop  open  surface  water  channel,  with  the  wall 
oscillating  at  a  frequency  of  T+  =  100  (2.67  Hz)  and  a  peak  to  peak  amplitude  of  2  cm.  Results  showed 
that  turbulent  velocity  fluctuations  were  considerably  reduced,  while  the  PIV  allowed  measurements  of  the 
physical  mechanisms  taking  place  in  the  boundary  layer.  The  main  result  of  these  experiments  was  that  the 
mean  velocity  increased  with  actuation,  and  was  accompanied  by  a  decrease  in  turbulence  production. 

Pang  and  Choi  [14]  recently  published  experiments  conducted  in  an  open  surface  water  channel  using 
Lorentz  force  actuation.  Actuators  cover  an  area  ~  300  mm  x  230  mm  where  a  high  conductivity  solution 
{CuSOa)  was  injected  at  the  surface  to  create  an  electrically  conductive  environment.  The  maximum  ve¬ 
locity  in  the  mean  profile  was  found  to  have  increased  approximately  66%  with  actuation,  when  u+  was 
nondimensionalized  using  its  own  friction  velocity  ( uT ). 

1.2  Review  of  scaling  arguments 

A  scaling  argument  was  developed  by  Breuer,  Park  and  Henoch  [1]  for  which  a  brief  review  will  be  presented 
here.  In  Lorentz  forcing,  in  order  for  the  control  to  be  effective,  the  force  must  be  able  to  induce  a  spanwise 
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velocity  in  the  fluid  of  approximately  uT.  The  Stuart  number,  St,  is  defined  as  the  ratio  of  the  Lorentz 
inertia,  JBa,  to  the  fluid  inertia,  pu\,  and  should  be  approximately  equal  to  1.  The  current  density,  J,  is 
defined  as  the  total  current  divided  by  the  total  wetted  area.  Thus  an  expression  can  be  obtained  for  the 
optimal  Stuart  number  at  which  Lorentz  forcing  will  be  effective: 


St  = 


JBa 
pul  ' 


(1.2) 


Defining  a  and  T  in  terms  of  turbulent  wall  units:  a  =  a+v/uT\T  =  T+v/ul  allows  us  to  redefine  the 
optimal  Stuart  number  as 


St  = 


a 

utT 


(1.3) 


We  can  extend  this  scaling  argument  slightly  further  for  the  current  experiments,  in  which  the  electrode 
spacing  a  has  a  fixed  size,  but  is  nevertheless  used  in  flows  with  different  Reynolds  numbers.  Thus  its 
“apparent”  size,  a+,  increases  as  the  Reynolds  number  of  the  flow  increases.  It  is  trivial  to  re-write  the 
definition  of  the  Stuart  number  (1.3)  as 


St  = 


a  Re* 

hY+ 


(1.4) 


where  h  is  the  channel  half-height  and  Re*  is  the  turbulent  Reynolds  number:  Re*  —  urhjv.  Thus  we 
see  that,  for  a  fixed  value  of  a,  we  predict  that  the  optimal  value  of  the  period  will  also  rise  linearly  with 
Reynolds  number,  and  that  for  a  fixed  non-dimensional  period,  the  optimal  Stuart  number  will  also  scale  with 
Re*.  This  gives  a  dimensionless  parameter  T+St/Re*  that  for  optimum  forcing  remains  constant,  regardless 
of  flow  rate.  The  experimental  data  confirms  these  predictions  both  for  the  drag  reduction  experiments 
performed  by  Breuer,  et  al.  [1]  as  well  as  the  flow  measurements  reported  here. 


1.3  Current  Experiments 

Experiments  reported  by  Breuer  et  al.  [l]  provided  a  basis  for  the  research  presented  here.  In  these  ex¬ 
periments,  Lorentz  force  actuators  were  assembled  together  to  form  a  large  plate  from  which  direct  drag 
measurements  were  obtained.  The  actuator  was  designed  for  robustness  in  the  harsh  saltwater  environment. 
The  electrodes  were  tested  extensively  for  effectiveness,  and  a  polymer  coating  was  found  to  decrease  both 
electrolytic  effects  during  actuation  as  well  as  corrosion  effects.  Spanwise  oscillation  as  well  as  travelling 
wave  actuation  were  tested.  Testing  was  somewhat  limited,  as  corrosion  and  degradation  of  the  actuators 
limited  experiment  duration.  Because  of  this,  design  changes  were  made  and  implemented  in  the  actuator 
fabrication  to  help  prevent  corrosion  and  increase  experimental  longevity.  The  experimental  procedure  and 
apparatus  was  also  changed  to  achieve  more  accurate  results  in  the  near  wall  region.  The  drag  measurements 
did  not  include  flow  data,  so  PIV  was  implemented  to  obtain  velocity  fields.  The  range  of  forcing  parameters 
was  expanded,  as  well  as  sampling  from  a  broader  range  of  Reynolds  numbers.  In  addition  to  providing  flow 
visualization,  the  data  further  confirms  the  scaling  arguments  presented  for  optimum  Lorentz  force  control. 
Experiments  were  run  to  study  the  effects  of  different  phases  of  the  forcing  cycle,  and  data  analysis  was 
expanded  to  include  integrated  actuation  data  for  both  low  and  high  Reynolds  numbers,  Ruu  and  Ruv  two 
point  correlations,  and  conditional  averaging. 
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Chapter  2 


Experimental  Setup  and  Procedure 


2.1  Actuator  Design  and  Fabrication 

The  actuator  design  is  based  on  the  actuator  developed  and  tested  by  Henoch  and  Stace  [15]  and  Breuer 
et  al.  [1]  which  used  arrays  of  magnets  and  electrodes  to  generate  a  uniform  Lorentz  force  acting  in  the 
streamwise  direction  (figure  2.1).  The  actuator  used  in  the  current  experiments  rotates  the  Henoch  and  Stace 
design  90  degrees,  so  the  Lorentz  force  now  acts  in  the  spanwise  direction,  perpendicular  to  the  core  flow. 
Linear  arrays  of  permanent  magnets  are  interlaced  with  a  commercially  fabricated  printed  circuit  board  to 
generate  the  necessary  fields. 

The  segmentation  of  the  electrodes  in  the  spanwise  direction  allows  multiple  forcing  regimens.  Spanwise 
oscillation,  simulated  by  Berger  et  al.  [16],  can  be  achieved  by  actuating  all  of  the  spanwise  columns  in  unison 
(i.e.  with  zero  phase  shift).  Operating  the  electrodes  with  the  same  amplitude  and  time-signal  but  including 
a  spanwise  phase  shift  generates  a  travelling  wave,  explored  computationally  by  Du  and  Karniadakis  [17]. 
They  found  that  actuation  via  a  transverse  travelling  wave  eliminated  near-wall  streaks,  reducing  shear  stress 
by  more  than  30%.  The  current  experiments  explore  spanwise  oscillation  only,  although  the  circuit  boards 
have  the  capacity  to  utilize  16  independent  channels,  making  travelling  wave  actuation  a  possibility. 

Based  on  the  experiments  run  by  Breuer  et  al.  [1]  some  design  modifications  have  been  made.  Previously, 
the  magnets  were  aligned  on  the  flux  plate  using  stainless  steel  pins.  The  installation  of  these  pins  was  both 
difficult  and  time-consuming.  To  improve  fabrication  time  as  well  as  reduce  corrosion,  plastic  dowel  pins  were 
substituted.  After  removal  of  the  actuators  from  the  previous  experiments,  it  was  found  that  the  tiles  were 
most  susceptible  to  saltwater  penetration  through  the  side  walls.  To  help  prevent  this  in  future  experiments, 
a  1/8”  epoxy  border  was  created  on  all  sides  of  the  actuator  (excluding  the  top  surface).  The  Nd-Bo-Fe 
permanent  magnets  used  were  also  found  to  have  suffered  corrosion  damage,  due  to  a  non-uniform  epoxy 
coating.  New  magnets  were  obtained,  which  were  both  nickel  plated  and  epoxy  coated  to  help  maintain  their 
integrity  throughout  the  course  of  the  experiments. 

Fabrication  of  the  actuators  was  a  time-consuming,  multi-step  process.  The  base  of  the  fixture  is  a 
12”  x  12”  x  1/8”  low  carbon  steel  plate.  The  plates  met  a  flatness  spec  of  0.008”  in  all  directions.  Low 
carbon  steel  was  chosen  for  its  machinability,  as  over  eight  hundred  holes  were  required  for  alignment  pins, 
epoxy  seepage,  and  structural  attachment.  Neodymium-Boron-Iron  (45  M-O)  permanent  magnets  (Adams 
Magnetics,  Elizabethtown  KY)  were  dual  coated  with  nickel  and  epoxy  to  help  prevent  degradation  in  the 
harsh  saltwater  environment.  The  magnets  have  a  1/8”  square  cross  section  and  a  length  of  2”.  Longer 
lengths  would  have  made  assembly  much  easier;  however,  the  strength  of  the  magnetic  field  as  well  as 
the  brittleness  of  the  magnets  were  limiting  factors.  The  magnets  were  aligned  in  rows  in  the  spanwise 
direction  using  plastic  dowel  pins.  These  3/32”  diameter,  1/8”  long  pins  were  inserted  individually  into  the 
plate  in  pre-drilled  holes,  making  the  assembly  very  tedious.  The  magnets  were  then  attached  to  the  plate 
and  maneuvered  into  place  against  the  plastic  pins,  giving  a  very  precise  array.  A  composite  photograph 
containing  both  the  magnet-pin  structure  as  well  as  the  electrode  layer  is  shown  in  figure  2.2. 

The  electrode  layer  consists  of  a  commercially  fabricated  two-layer  printed  circuit  board  (Multi-Lab,  Inc. 
Manchester,  MA).  The  electrodes  protrude  slightly  above  the  surface,  but  are  hydrodynamically  smooth  (step 
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Figure  2.2:  Composite  photograph  of  Lore: 
pins  and  magnets)  and  the  electrode  layer  [; 
force  is  from  top  to  bottom.  The 


Figure  2.1:  Schematic  of  the  Lorentz  Force 
force  acts  in  the  spanwise  direction,  out  of 
view  [1]. 


Lorentz  force 
directon 


Electrode  pair 


Figure  2.3:  Photographs  showing  (a)  PCB  in  aluminum  mold  with  tape  border  and  1/8”  spacer  inserted  into 
corner  and  (b)  finished  plate  with  epoxy  coating  and  “levelling  feet”  embedded  on  the  back  surface.  Holes 
drilled  through  the  flux  plate  to  allow  seepage  are  also  visible. 
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Figure  2.4:  Schematic  of  flow  facility,  showing  the  back  pressure  and  water  supply  tanks,  development 
section,  test  area  with  actuator  and  return  flow  path. 

less  than  If).  To  increase  the  effective  electrode  separation  the  electrodes  are  designed  with  a  “cut  back” 
between  spanwise  columns.  The  surface  of  the  circuit  board  was  screen  printed  with  a  polymer  ink  to  both 
reduce  the  electrolytic  effects  during  actuation  and  prevent  corrosion.  Thirty-two  wires  were  soldered  to  each 
board  with  a  length  of  21  feet  so  that  each  wire  could  be  directly  connected  to  the  amplifiers,  alleviating  the 
need  for  cumbersome  bus  lines. 

A  mold  was  designed  (figure  2.3)  to  facilitate  the  assembly  of  the  magnet/flux  plate  fixture  to  the  circuit 
board.  The  inner  dimensions  of  the  Aluminum  mold  were  machined  to  be  121/4”  square,  to  allow  for  the 
epoxy  border  around  the  outside.  The  walls  of  the  mold  are  detachable,  giving  easy  removal  after  the  epoxy 
has  cured.  To  maintain  a  constant  mold  size,  a  plastic  template  precision  machined  to  the  correct  dimensions 
was  used  to  size  the  mold  before  each  use.  All  of  the  inner  surfaces  were  then  coated  with  a  silicone  spray,  to 
prevent  sticking  and  subsequent  cracking  of  the  epoxy  upon  removal.  The  top  (actuated)  surface  of  the  PCB 
was  covered  with  a  sheet  of  plastic  and  masking  tape  to  prevent  seepage  of  epoxy  onto  the  electrodes.  A  1/8” 
border  of  masking  tape  was  allowed  on  the  outside  of  the  circuit  board  to  center  it  in  the  mold.  Levelling 
feet  were  attached  to  the  bottom  surface  of  the  steel  plate  using  brass  set  screws.  These  1”  diameter,  1/2” 
thick  polypropylene  discs  were  used  to  assure  flatness  of  the  plates  when  installed  into  the  channel,  as  well 
as  aiding  in  the  removal  of  the  plates  from  the  mold.  The  magnet /flux  plate  fixture  was  suspended  above 
the  mold  containing  the  circuit  board  while  the  wires  were  individually  fed  through  the  assigned  holes. 

The  PCB  was  attached  to  the  magnet/flux  plate  fixture  using  a  commercially  available  epoxy  (Emerson 
and  Cuming  Stycast  2057  black  epoxy),  chosen  because  it  is  waterproof  and  corrosion  resistant.  Once  the 
circuit  board  was  mounted  in  the  mold  and  the  wires  were  fed  through  the  steel  base  plate,  the  epoxy  was 
mixed  continuously  with  the  hardener  for  two  minutes  to  ensure  continuity.  The  mixture  was  poured  over 
the  PCB,  and  the  flux  plate  was  lowered  in  increments  to  make  sure  none  of  the  wires  were  crimped  and 
subsequently  damaged  when  the  actuator  was  clamped.  To  make  sure  the  actuator  plate  was  centered  in 
the  mold,  four  90  degree  1/8”  L-bracket  spacers  were  placed  in  the  corners.  The  spacers  were  removed  only 
after  the  plate  was  completely  lowered  into  the  mold  and  clamps  applied  holding  the  levelling  feet  firmly  to 
the  mold.  The  remainder  of  the  epoxy  mixture  was  poured  and  evened,  and  was  given  a  full  24  hours  to 
cure.  After  complete  hardening  the  mold  was  carefully  removed  piece  by  piece  to  avoid  cracking  the  plate, 
and  all  actuators  were  checked  for  continuity. 
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Figure  2.5:  Schematic  showing  side  and  top  views  of  PIV  setup.  A  laser  sheet  in  the  x-y  plane  is  produced 
by  an  Nd:YAG  laser  mounted  above  the  channel.  The  image  size  is  approximately  9  x  11mm. 


Parameter 

Flow  Rate  1 

Flow  Rate  2 

Flow  Rate  3 

Flow  Rate  4 

uT  [m/s] 
t*  [ms] 
l*  [/i  m] 

Re* 

Frequencies  [Hz] 
Amplitudes  [Vo_p] 
St 

T+ 

0.0064 

26.0 

166.10 

144 

(0.3  -  3.2) 
(0.5  -  3.5) 
(2.9  -19.0) 
(12  -  123) 

0.0132 

5.7 

74.9 

305 

(0.3  -  1.6) 
(3.0  -  5.0) 
(3.5  -  6.8) 
(104  -  564) 

0.0147 

4.5 

67.5 

330 

(0.8  -  4.1) 
(2.0  -  5.5) 
(2.2  -  5.7) 
(52  -  258) 

Table  2.1:  Forcing  frequencies,  amplitudes,  dimensionless  numbers  and  flow  parameters  for  the  four  flow 
rates  used. 

2.2  Flow  Facility 

The  flow  facility  (figure  2.4)  at  Brown  University  is  a  constant  pressure  channel  that  can  support  both 
saltwater  and  freshwater  experiments,  with  flow  rates  between  Re *  =  100  and  Re *  =  800.  The  entire 
channel  is  constructed  of  clear  acrylic  for  ease  of  PIV  measurement  as  well  as  corrosion  resistance  and  is 
built  on  80-20  supports.  Flow  is  pumped  from  the  collection  tank  through  a  series  of  four  pumps  up  to  the 
accumulation  tank.  Flow  then  continues  into  the  development  section  through  a  honeycomb  and  over  a  strip 
of  velcro  that  forces  early  transition  to  turbulence.  The  development  section  is  3m  long,  lm  wide,  and  24mm 
high.  The  flow  then  enters  the  specially  designed  test  section. 

The  floor  of  the  test  section  is  a  removable  table  that  allows  installation  and  adjustment  of  multiple 
actuators.  Snorkels  in  the  roof  of  the  channel  as  well  as  a  “fence”  along  the  side  of  the  test  section  permit 
the  wires  to  be  fed  out  of  the  channel  to  the  amplifiers  without  disturbing  the  main  flow.  Portholes  in  the 
roof  give  easy  access  for  calibration  and  cleaning.  Salt  was  added  to  a  mass  concentration  of  3.5%  giving  a 
conductivity  of  4.5  Siemens/m,  closely  matching  that  of  seawater.  The  channel  was  then  run  for  eight  hours 
to  provide  adequate  time  for  mixing. 

Numerous  control  and  monitoring  devices  have  been  installed  for  precise  control  and  repeatability  of 
experiments.  It  was  found  that  the  temperature  increase  due  to  work  added  to  the  system  by  the  pumps 
was  approximately  three  degrees  per  hour,  significantly  changing  the  viscosity.  A  heat  exchanger  (NESLAB 
HX-150DD)  was  used  in  the  collection  tank  to  reduce  the  temperature  variation  to  within  0.2  degrees  over  a 
three  hour  experiment.  A  valve  system  allows  flow  to  bypass  the  main  channel,  giving  more  precise  control 
of  flow  rate.  The  flow  rate  was  continually  monitored  during  testing  using  a  fixed-depth  turbine  flow  sensor 
(Seametrics). 

Trapped  air  under  the  plates  caused  bubbles  during  PIV  measurements,  decreasing  the  clarity  of  the 
image.  The  surface  of  the  actuator  tiles  was  cleaned  before  each  experiment  and  a  specially  designed  wiper 
allowed  for  “sweeping”  of  the  test  surface  as  well  as  the  roof  of  the  channel  during  testing.  The  temperature, 
flow  rate  and  pressure  head  were  monitored  and  recorded  during  each  actuation.  Sinusoidal  waveforms  were 
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Figure  2.6:  Experimentally  measured  voltage  divided  by  current  plotted  with  actuation  frequency  (T+).  V/I 
remains  approximately  constant  for  all  experiments. 


produced  using  a  Stanford  Research  Systems  Waveform  Generator  and  a  16-channel  custom  built  amplifier. 
Power  is  supplied  by  2-10V  power  supplies  (Lambda),  and  current  to  the  actuators  was  continuously  mon¬ 
itored  using  an  open  loop  current  sensor  (Sypris,  IHA-25).  Figure  2.6  shows  the  experimentally  measured 
voltage  divided  by  current  plotted  with  the  forcing  period,  T+.  As  expected,  the  value  remains  constant  for 
all  experiments. 


2.3  Particle  Image  Velocimetry  System 

A  Particle  Image  Velocimetry  (PIV)  System  was  used  to  obtain  velocity  data.  Pairs  of  images  were  acquired 
from  which  full  velocity  fields  were  calculated.  A  schematic  showing  the  PIV  setup  is  presented  in  figure  2.5. 
A  15  mW  Nd:YAG  laser  is  mounted  on  a  traverse  which  moves  in  the  spanwise  direction,  with  attached  optics 
that  create  a  laser  sheet  in  the  x  —  y  plane.  The  laser  and  the  camera  can  be  moved  such  that  measurements 
can  also  be  taken  in  the  x  —  z  plane.  A  3D  precision  stage  with  micrometer  adjustment  capabilities  holds 
the  camera  and  allows  for  fine  corrections  during  testing.  An  Infinity  lens  (InFocus  Model  KC)  which  has 
both  a  large  working  distance  (ss  350mm)  as  well  as  a  large  field  of  view  allowed  magnified  images  near  the 
wall.  PIV  software  (IDT  pro  VISION,  release  2.01)  was  used  to  process  the  images  after  acquisition.  Image 
calibration  was  achieved  using  a  1”  Aluminum  block  drilled  with  holes  in  a  2mm  square  grid.  The  laser 
sheet  was  located  approximately  2cm  upstream  from  the  end  of  the  actuator  plates,  and  was  centered  on  the 
actuator  in  the  spanwise  direction. 

2.4  Experimental  Procedure 

Before  each  data  acquisition,  the  setup  was  allowed  to  run  for  approximately  30  minutes.  This  allowed 
the  flow  rate  to  become  constant  as  well  as  any  bubbles  to  be  flushed  out  of  the  channel.  Bubbles  on  the 
actuated  surface  scattered  the  laser  light  and  blurred  the  image,  and  proved  difficult  to  get  rid  of.  The  roof 
of  the  channel  also  gathered  bubbles  easily,  distorting  the  laser  sheet.  A  specially  designed  ‘bubble  scraper’ 
allowed  the  surfaces  to  be  cleaned  during  testing.  This  allowed  the  integrity  of  the  image  to  be  maintained 
through  the  course  of  the  experiment,  without  altering  the  flow  rate  or  pressure  head. 

After  the  system  was  allowed  to  come  to  equilibrium,  an  initial  baseline  case  was  taken.  The  data  was 
immediately  processed  using  the  PIV  software  to  assure  accurate  results.  Barring  any  abnormalities  in  the 
flow,  the  test  matrix  was  then  run.  In  between  each  case,  the  channel  was  given  30  seconds  to  allow  the 
new  forcing  scheme  to  propagate  through  the  entire  length  of  the  test  section.  Baseline  cases  were  taken 
throughout  the  course  of  the  experiment  to  ensure  no  drifting  occurred. 
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Figure  2.7:  PIV  baseline  channel  velocity  profiles  showing  both  mean  flow  and  fluctuations.  Profiles  demon¬ 
strate  fully  developed  turbulent  flow  at  Re*  =  167.  All  quantities  are  normalized  by  the  friction  velocity, 


In  early  experiments,  glare  from  the  electrodes  made  obtaining  accurate  near  wall  data  difficult.  By 
adjusting  the  contrast  and  laser  intensity,  consistent  measurements  were  made  to  y+  ~  4.  A  smaller  mesh 
size  and  corresponding  interrogation  area  were  used  close  to  the  wall,  allowing  a  more  accurate  correlation. 
The  1360x1024  pixel  image  was  covered  with  a  mesh  that  created  a  velocity  field  with  approximately  50x75 
vectors.  The  camera  was  rotated  90  degrees  to  allow  the  longer  dimension  to  be  perpendicular  to  the  wall. 
The  typical  PIV  interrogation  area  used  was  36  pixels.  Time  separation  between  the  laser  pulses  for  the 
lower  Reynolds  number  was  300/i.s,  compared  to  150/x.s  for  the  higher  velocity  flow.  Each  point  in  the  velocity 
profile  was  obtained  by  averaging  over  approximately  50  vectors  in  the  x-direction  and  200  ensemble  images. 
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Chapter  3 


Results  and  Discussion 


3.1  Mean  Flow 

A  series  of  experiments  were  run  at  two  different  Reynolds  numbers.  Figure  2.7  shows  baseline  channel 
velocity  profiles  for  the  lower  Reynolds  number,  Re*  =  167.  The  mean  flow  is  on  the  left  with  the  fluctuations, 
u',  v',  and  —u'v'  on  the  right,  showing  a  well  behaved  and  fully  developed  turbulent  flow.  With  the  previously 
mentioned  adjustments  to  the  experimental  procedure,  accurate  measurements  were  made  to  within  y+  =  4. 
A  second  set  of  experiments  was  run  at  a  higher  velocity  to  study  the  effects  of  forcing  with  increased 
Reynolds  number.  Baseline  profiles  are  presented  in  figure  2.8.  All  baseline  as  well  as  actuated  cases  have 
been  normalized  by  the  friction  velocity  ( uT )  of  the  unactuated  case.  Since  the  water  channel  is  driven  with 
a  constant  pressure  head,  uT  remains  constant  for  all  cases.  Initially,  uT  was  calculated  using  u+  —  y+ 
in  the  near  wall  region,  and  this  gave  a  rough  approximation.  By  matching  the  unactuated  mean  profiles 
with  turbulence  statistics  from  a  Direct  Numerical  Simulation  (DNS)  [2],  uT  was  re-scaled  to  give  a  more 
accurate  fit.  The  experimental  data  scales  well  with  the  numerical  calculations,  and  the  results  are  presented 
in  figure  3.1.  The  close  fit  not  only  shows  good  agreement  between  numerical  simulations  and  experimental 
results,  but  also  provides  a  very  accurate  method  for  determining  uT. 


3.2  First  Order  Statistics 

3.2.1  Actuation  Data 

Figures  3.2  through  3.5  show  profiles  of  the  mean  velocity,  Reynolds  stress,  and  turbulence  production 
terms.  For  each  figure,  either  the  frequency  or  amplitude  is  held  constant  and  results  are  shown  for  a  sweep 
through  the  other  parameter.  In  all  figures  presented  here,  blue  circles  represent  the  baseline  (unactuated) 
case.  Figure  3.2  shows  a  frequency  sweep  for  a  fixed  amplitude  (4.5  V,  St  =  12.9)  at  the  lower  speed 
flow  (Re*  =  167).  The  red  circles  show  the  lowest  frequency,  green  circles  the  highest,  and  the  black  circles 
represent  the  frequency  that  was  found  to  produce  a  maximum  in  local  flow  rate  accompanied  by  a  minimum 
in  Reynolds  stress.  Actuated  cases  show  both  an  increase  in  flow  rate  as  well  as  a  decrease  in  Reynolds  stress, 
demonstrating  that  the  Lorentz  forcing  is  indeed  affecting  the  near  wall  flow. 

Given  the  profiles  for  the  mean  flow  and  Reynolds  stress,  we  can  calculate  the  turbulence  production 
term.  Taken  from  the  fluctuation  energy  equation,  the  turbulence  production  is  defined  as 

P  «  -W  •  dU/dy  (3.1) 

Figure  3.3  shows  the  Reynolds  stress  and  turbulence  production  profiles  for  the  cases  shown  above.  As 
expected,  a  decrease  in  production  accompanies  an  increase  in  the  mean  flow  and  decrease  in  Reynolds  stress. 
Production  is  calculated  to  y+  ~  100  as  in  the  mean  flow  and  Reynolds  stress  cases,  but  is  plotted  here  only 
to  y+  =  50  to  magnify  the  near  wall  region.  A  ‘sweet  spot’  at  T+St/Re*  =  2.9  was  found,  showing  the 
greatest  decrease  in  turbulence  production  accompanied  by  the  largest  increase  in  flow  rate. 
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Figure  3.1:  Baseline  unactuated  mean  profiles  plotted  with  numerical  turbulence  statistics  for  a  range  of 
Reynolds  numbers  from  Moser  et  al.  [2]. 


/  y* 


Figure  3.2:  Mean  flow  profiles  for  a  fixed  forcing  amplitude  with  changing  frequency  at  Re *  =  167.  The  fixed 
amplitude  that  was  chosen  was  the  highest  possible  given  the  actuator  design.  Actuation  both  increases  and 
decreases  the  fluctuations,  with  a  local  minimum  around  T+St/Re*  =  2.9. 
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Figure  3.3:  Changes  in  Reynolds  stress  and  turbulence  production  for  a  constant  amplitude  and  changing 
frequencies.  The  lowest  production  accompanies  the  highest  local  flow  rate  as  well  as  the  greatest  reduction 
in  shear  stress. 


Figure  3.4:  Effects  of  forcing  on  the  mean  flow  are  shown  for  a  higher  flow  rate  (Re*  =  330).  One  forcing 
frequency  and  three  amplitudes  are  plotted  withthe  baseline. 
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Figure  3.5:  Changes  in  Reynolds  stress  and  turbulence  production  at  Re*  =  330.  The  turbulence  production 
is  greatly  suppressed  around  the  ‘sweet  spot’  where  T+St/Re*  =  2.2. 


Figure  3.6:  Percent  change  in  flow  rate  and  production  terms  over  a  wide  range  of  forcing  frequencies  and 
amplitudes  for  a  flow  rate  of  Re*  =  167.  The  maximum  increase  in  flow  rate  is  accompanied  by  the  greatest 
decrease  in  turbulence  production,  near  the  ‘sweet  spot’. 


Figures  3.4  and  3.5  show  the  effects  of  forcing  for  a  higher  Reynolds  number  (Re*  =  330).  The  actuated 
cases  cause  a  change  in  the  flow;  however,  the  magnitude  is  much  lower.  This  can  be  predicted  from  the 
scaling  arguments  made  earlier  by  Breuer,  et  al.  [1].  The  scaling  predicts  that  for  a  constant  optimum 
value  of  T+  St/Re* ,  when  the  Reynolds  number  is  roughly  doubled  the  optimum  forcing  comes  when  the 
quantity  T+St  is  also  doubled.  Due  to  experimental  restrictions  on  the  amount  of  current  that  can  be  run 
through  the  system,  it  was  not  possible  to  effectively  double  this  parameter.  This  being  said,  there  is  still  a 
significant  change  in  the  Reynolds  stress  as  well  as  the  production  term.  The  greatest  decrease  in  turbulence 
production  comes  at  the  largest  forcing  amplitude  which  correlates  to  a  T+St/Re*  ~  2.2.  Since  a  minimum 
in  production  has  not  yet  been  reached,  it  is  foreseeable  that  increasing  the  current  through  the  system  can 
further  decrease  the  turbulence  production  and  increase  the  mean  flow. 

3.2.2  Integrated  Flow  Rate  and  Production 

The  actuated  results  for  each  Reynolds  number  are  effectively  summarized  in  figures  3.6  and  3.7.  For  each 
actuated  case,  the  mean  flow  and  production  profiles  were  integrated  to  y+  ~  75  to  find  a  local  flow  rate  and 
integrated  production  value.  The  percent  change  from  the  baseline  for  the  flow  rate  is  shown  on  the  left, 


13 


Figure  3.7:  Percent  change  in  flow  rate  and  production  terms  over  a  wide  range  of  forcing  frequencies  and 
amplitudes  for  a  flow  rate  of  Re *  =  334.  The  data  shown  represents  five  forcing  frequencies  and  five  forcing 
amplitudes. 


with  the  accompanying  change  in  production  on  the  right.  All  forcing  amplitudes  are  presented  across  the 
horizontal  axis,  spanning  St  —  6.4  —  17.9.  Forcing  frequencies  are  represented  on  the  vertical  axis,  covering 
T+  =  13  —  265.  A  maximum  increase  in  flow  rate,  coupled  with  a  minimum  in  turbulence  production  can  be 
seen  at  T+  =  32  and  a  St  =  17.  This  correlates  to  a  T+ St/Re*  =  2.9,  which  agrees  well  with  the  optimum 
conditions  suggested  by  the  mean  flow  and  Reynolds  stress  plots  presented  earlier.  It  is  also  observed  that 
there  are  regions  where  the  flow  rate  decreases  and  the  turbulence  production  increases,  which  shows  that 
Lorentz  forcing  does  indeed  have  a  significant  impact  on  the  flow. 

Integrated  mean  flow  and  turbulence  production  data  for  a  Re*  —  330  is  presented  in  figure  3.7.  The 
flow  rate  increases  for  all  actuated  cases,  which  is  expected  considering  there  were  no  forcing  parameters  that 
decreased  the  mean  flow  or  increased  the  Reynolds  stress.  This  suggests  that  optimum  testing  conditions 
have  not  been  reached.  As  stated  previously,  the  actuators  limit  the  amount  of  current  that  can  be  run 
through  the  system.  Although  all  actuated  cases  decrease  the  production,  the  maximum  change  is  still  20% 
less  than  for  the  lower  Reynolds  number  flow,  suggesting  higher  forcing  amplitudes  would  be  more  effective. 

Figure  3.8  represents  all  of  the  data  taken  in  the  current  experiments,  scaled  by  the  dimensionless 
parameter,  T+St/Re*.  The  percent  change  from  the  unactuated  case  for  both  flow  rate  and  production 
are  shown.  Results  from  four  separate  experiments  covering  two  ranges  of  Reynolds  numbers  are  presented. 
Both  the  flow  rate  and  turbulence  production  data  scale  well  using  the  dimensionless  ratio,  and  suggest 
an  optimum  forcing  condition  of  T+St/Re*  «  2  —  3.  Lorentz  forcing  has  less  of  an  impact  on  the  higher 
Reynolds  number  flow,  showing  a  smaller  change  in  flow  rate  and  production.  A  discrepancy  exists  between 
the  values  of  T+  and  St  numbers  presented  here,  and  those  which  appear  in  the  Breuer  et  al.  [1]  results. 
The  difference  is  most  likely  a  result  of  the  calculation  of  the  current  density,  J.  In  the  current  experiments, 
current  was  rigorously  monitored  for  each  case  and  the  recorded  number  used  directly  in  calculations.  The 
current  density  in  the  previous  experiments  was  most  likely  derived  from  the  forcing  amplitude,  and  thus 
was  not  as  accurate. 

3.2.3  Phase  Synchronization 

A  separate  experiment  was  run  to  determine  the  effects,  if  any,  of  the  different  phases  of  the  forcing  cycle  on 
the  mean  flow  and  turbulence  production.  As  mentioned  in  the  experimental  procedure,  a  sinusoidal  wave 
was  generated  and  amplified  to  drive  the  actuators.  Data  was  taken  at  90°,  180°,  270°,  and  360°  and  plotted 
in  comparison  to  the  baseline.  Figure  3.9  shows  the  mean  profiles  and  the  Reynolds  stress  for  the  actuated 
case  chosen  ( T+St/Re *  =  4.7).  At  the  time  the  in-phase  measurements  were  taken,  the  ‘sweet  spot’  of 
actuation  had  not  yet  been  determined,  and  thus  a  high  voltage,  low  frequency  case  was  chosen  at  random. 
The  Reynolds  stress  and  therefore  the  production  for  some  phases  go  negative,  which  was  not  observed  in 
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Figure  3.8:  Changes  in  the  turbulence  production  and  flow  rate  as  functions  of  the  scaling  parameter, 
StT+  /Re*.  Data  from  four  different  Reynolds  numbers  are  shown  and  a  wide  range  of  forcing  amplitudes 
and  frequencies.  Data  collapses  well,  showing  a  minimum  in  production  and  maximum  in  flow  rate  at 
T+St/Re*  between  2  and  3. 


Figure  3.9:  Mean  and  Reynolds  Shear  Stress  taken  at  four  different  phases  of  the  forcing  cycle  for  a  Re*  =  200. 
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Figure  3.10:  Turbulence  production  calculated  for  different  phases  of  the  forcing  cycle. 


Figure  3.11:  RuniyuUz)  and  RVv(yi,  1/2)  at  the  location  of  maximum  turbulence  production,  y+  =  12.  In  all 
actuated  cases,  the  correlation  decreases,  with  a  sweet  spot  at  approximately  StT+  /Re*  =  2.9 


any  previous  data.  The  variation  between  the  cases  is  similar  to  that  seen  between  baseline  cases  through 
the  course  of  one  experiment.  The  source  of  the  negative  values  could  possibly  be  attributed  to  experimental 
error;  the  experiments  are  inconclusive  and  will  be  run  again  to  find  a  definitive  explanation. 


3.3  Second  Order  Statistics 

3.3.1  Two-point  correlations 

Correlation  coefficients  were  computed  for  each  of  actuated  cases  run.  The  correlation  coefficient  is  defined 
as 


=  (3.2) 

v°2(^)  •  v 62  (2/2) 

Ruu  5  Rvv  j  and  Ruv  were  calculated  for  no  control  as  well  as  all  actuated  cases  at  specific  values  of 
y+ .  Figure  3.11  shows  Ruu  and  Rvv  at  the  location  of  maximum  turbulence  production  for  the  lower 
velocity  flow  (Re*  =  167).  All  actuated  cases  show  a  lower  correlation  than  the  baseline  case  (blue  circles). 
The  greatest  change  comes  in  the  near  wall  region,  due  to  interaction  with  coherent  structures.  In  this 
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Figure  3.12:  RUu{yi-,yi)  in  the  log  layer  for  (a)  Re*  =  167  at  y+  —  76  and  (b)  Re*  =  330  at  y+  =  122. 


region,  the  T+St/Re*  that  produced  the  greatest  increase  in  flow  rate  and  decrease  in  turbulence  production 
(T+ St/Re*  —  2.9)  also  represents  the  greatest  decrease  in  correlation.  In  the  Rvv  correlation,  the  actuated 
cases  also  show  a  deviation  from  the  baseline.  The  lowest  frequency  (red  circles,  T+St/Re*  =  6.7)  shows  the 
least  impact  on  the  flow,  while  the  higher  frequencies  show  approximately  the  same  reduction  in  correlation. 
A  lower  correlation  of  Ruu  implies  either  smaller  coherent  structures,  fewer  structures,  or  less  coherence  in 
the  structure.  It  is  inferred  that  the  near  wall  disturbances  created  by  the  Lorentz  forcing  cause  the  lowest 
correlation  to  occur  when  y+  <  10,  as  it  is  in  this  region  that  the  forcing  has  the  greatest  impact.  Figure  3.12 
shows  the  correlations  calculated  for  locations  in  the  log  layer  for  both  the  low  and  high  Reynolds  number 
flows.  The  effect  of  the  forcing  on  the  correlation  decreases  with  distance  from  the  wall  as  expected,  given 
the  penetration  depth,  A+  «  9  for  the  lower  Reynolds  flow. 

3.3.2  Conditional  Averaging 

Based  on  the  work  of  Johansson  et  al.  [18]  probabilities  for  the  Reynolds  shear  stress  were  calculated.  For 
each  acquisition,  the  u'  and  v'  were  determined  by  subtracting  the  mean  flow  from  the  non-averaged  total 
flow.  Using  these  values,  an  instantaneous  value  of  —  u'v'  was  determined.  Dividing  this  instantaneous 
value  by  the  averaged  value,  —  u'v' ,  gave  an  indicator  of  turbulent  events.  This  probability,  P,  is  plotted 
against  the  threshold  value  k  in  figure  3.13  for  both  high  and  low  Reynolds  flow.  The  baseline  case  is  again 
represented  with  blue  circles,  with  the  other  colors  representing  a  frequency  sweep  at  a  constant  amplitude. 
For  the  low  Reynolds  flow  (left  figure)  the  number  of  high  intensity  events  decreases  at  optimum  forcing 
(black  circles),  as  expected.  The  high  frequency  case  (red  circles)  which  previously  showed  a  decrease  in  flow 
rate  and  an  increase  in  turbulence  production,  here  demonstrates  an  increase  in  turbulence  intensity  events. 
In  the  case  of  the  high  Reynolds  number  flow,  the  actuated  cases  show  no  difference  from  the  baseline  for 
a  threshold  parameter  k  less  than  three.  For  values  greater  than  three  all  actuated  cases  produce  a  greater 
number  of  high  intensity  events,  which  was  not  the  expected  result.  It  is  proposed  that  with  a  greater 
St  number  (forcing  voltage),  a  greater  difference  in  the  number  of  high  intensity  events  will  be  observed. 
Further  analysis  will  be  conducted  for  a  greater  understanding  of  these  results. 
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Figure  3.13:  Probabilities  plotted  as  a  function  of  the  threshold  parameter,  k,  for  a  frequency  sweep  at  a 
constant  forcing  amplitude  for  both  low  and  high  Reynolds  number  flow. 
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Chapter  4 

Conclusions 


Lorentz  force  actuators  were  tested  in  low  Reynolds  number  saltwater  channel  flow  to  study  the  effects  of  near 
wall  fluid  perturbations.  Robust  actuators  were  designed,  fabricated  and  tested  and  mounted  in  the  Brown 
University  flow  facility.  Particle  Image  Velocimetry  was  used  to  produce  flow  profiles  for  the  mean  flow  as 
well  as  fluctuations  for  baseline  and  actuated  cases.  A  range  of  forcing  frequencies  and  amplitudes  were  tested 
based  on  the  optimum  forcing  predicted  by  the  dimensionless  parameter  T+ St/Re*  through  two  different 
Reynolds  number  ranges.  From  the  mean  flow  profiles,  a  local  flow  rate  was  determined  and  found  to  increase 
greatly  with  the  optimal  forcing.  Similarly,  for  optimum  forcing  parameters,  the  turbulence  production  was 
calculated  and  was  found  to  decrease  significantly  with  Lorentz  force  control.  Control  was  greater  with  the 
lower  Reynolds  number  flow,  as  predicted,  as  experimental  constraints  prevented  scaling  the  T+  and  St 
numbers  properly  with  the  increase  in  Re*  number.  Both  the  local  flow  rate  and  turbulence  production 
were  integrated  over  the  entire  range  of  scaling  parameters,  and  broad  areas  of  both  positive  and  negative 
control  were  seen.  The  percent  change  from  the  baseline  for  both  flow  rate  and  production  was  plotted  as  a 
function  of  the  dimensionless  parameter  T+ St/Re*,  and  the  optimal  forcing  was  seen  at  T+ St/Re*  «  2  —  3. 
The  experimental  data  taken  in  phase  with  the  sinusoidal  forcing  proved  to  be  inconclusive,  and  further 
experiments  are  needed  to  draw  any  definitive  conclusions.  Two  point  correlations  were  computed  for  all  of 
the  actuated  cases,  inferring  a  reduction  in  the  size  of  the  near  wall  coherent  structures  with  Lorentz  forcing. 
Conditional  sampling  was  performed  and  low  Reynolds  number  results  were  in  agreement  with  the  actuated 
data  and  correlation  findings.  Future  research  will  include  taking  data  in  the  x  —  z  plane  to  study  the  near 
wall  structures,  repeating  the  phase-locked  forcing  data,  and  further  analyzing  the  conditional  averaging 
data.  Separation  control  using  the  Lorentz  force  has  been  proposed,  and  will  hopefully  be  investigated  in 
the  near  future. 
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